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Abstract 

We investigate the robustness of parity- and time-reversal (VT) symmetric phase in an TV-site 
lattice with position-dependent, parity-symmetric hopping function and a pair of imaginary, VT- 
symmetric impurities. We find that the "fragile" 'PT-symmetric phase in these lattices is stronger 
than its counterpart in a lattice with constant hopping. With an open system in mind, we explore 
the degrees of broken PT-symmetry and their signatures in single-particle wavepacket evolution. 
We predict that when the 'PT-symmetric impurities are closest to each other, the time evolution of 
a wavepacket in an even-iV lattice is remarkably different from that in an odd- N lattice. Our results 
suggest that PT-symmetry breaking in such lattices is accompanied by rich, hitherto unanticipated, 
phenomena. 



Introduction: Since their discovery by Bender and coworkers over a decade ago [U |2], non- 
Hermitian Hamiltonians Hp? that are symmetric with respect to combined parity (V) and 
time- reversal (T) operations have been extensively explored [31 H] . Such continuum Hamilto- 
nians usually feature a non-Hermitian, PT-symmetric potential, V(x) = V*(— x) ^ V*(x). 
Although not Hermitian, Hpt 7^ H PT , they have purely real eigenvalues over a range of 
parameters. Traditionally, the region of parameters where the spectrum of the Hamiltonian 
is real and, therefore, its eigenfunctions are simultaneous eigenfunctions of the combined 
PT-operation, is defined as the PT-symmetric region. When an eigenvalue of the PT- 
symmetric Hamiltonian becomes complex, E ^ E*, since time- reversal is an anti-linear 
operation, it follows that the corresponding eigenfunction is not invariant under the PT- 
operation, [Hpt,VT]iPe{x) 7^ 0. Therefore, the number of eigenfunctions that break the 
PT-symmetry is equal to the number of complex eigenvalues. 

In the past few years, PT-symmetric lattice models have become a focal point of re- 
search [IHTj. Lattice models, in general, are popular in physics due to their versatility [8], 
availability of exact solutions [9], the absence of an ultraviolet divergence [10], and the 
ability to capture counterintuitive phenomena that have no counterparts in the continuum 
theories [TT] . Recently, evanescently coupled waveguides [12] have emerged as a promising 
candidate for the realization of an ideal, one-dimensional lattice with tunable hopping [13J, 
disorder [13] , and non-Hermitian, on-site, impurity potentials [ToTHTj . The hopping between 
adjacent waveguides is controlled by the width of the barrier between the two, and the com- 
plex, on-site potential is proportional to the local refractive index, the imaginary part of 
which provides the loss or gain in the respective waveguide. The emergence of complex 
eigenvalues in such a pair of waveguides, via the attendant non-reciprocal behavior of light 
intensity, has been experimentally explored [TBI [T7]. A salient feature of a PT-symmetric 
lattice with (approximately) constant hopping is that its PT-symmetric phase is fragile ex- 
cept when the non-Hermitian, PT-symmetric impurities are closest or the farthest [51I71IT8]. 
The PT-symmetric phase becomes robust if, instead, a non-Hermitian, PT-symmetric hop- 
ping function is chosen [19] . However, the degrees of PT-symmetry breaking, defined as the 
fraction of eigenvalues that become complex, and their consequences to the properties of the 
system, have not been explored [2"U] . 

In this paper, we investigate the robustness of "PT-symmetric phase, and the degrees of 
PT-symmetry breaking in an iV-site lattice with a hopping function t a (k) = to[k(N — k)] a / 2 
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and a pair of PT-symmetric impurities ±27 at positions (mo, mo) respectively where mo = 
N + 1 — m is the mirror-symmetric counterpart of site mo- We define the PT-symmetric 
phase as robust (fragile) provided the critical impurity strength ^pt{N), below which all 
eigenvalues are real, is nonzero (zero) in the limit N — > oo. With coupled optical waveguides 
in mind [TBI ITT] , we explore the signatures of PT-symmetry breaking in the evolution of 
single-particle properties across the PT-symmetric phase boundary. 

Our main results are as follows: i) the PT-symmetric phase is robust when a > and 
the impurities are closest to each other, ii) in the fragile case, the scaling of 7p T (iV) depends 
on the fractional impurity position \i = tjiq/N and a, iii) in the robust case, single-particle 
wavepacket evolution across the PT-symmetric phase boundary is qualitatively different for 
even and odd N. Our results show that, based on the position of impurities and the form 
of the hopping function, signatures of the PT-symmetry breaking in lattices with N ^> 1 
are dramatically different from those in the experimentally explored N = 2 case [161 H7J- 
Tight-binding Model: We start with the Hamiltonian for a one-dimensional lattice 

N-l 

H a = - ^2 *«( ? ") ( a l+i a * + a \ a i+i) + h ( a L am o _ a k a mo) > C 1 ) 
i=i 

where a\(ak) is the creation (annihilation) operator for a state localized at site k, and 
1 < mo < N/2 is the position of the first impurity. We recall that for a constant hopping 
model, as 7/to is increased, the maximum number of complex eigenvalues is given by 2mo = 
2/iiV < N. This process is always sequential: the number of complex eigenvalues increases 
by two or four as 7 increases; the sole exception is the case with nearest neighbor impurities, 
m = N/2, where all eigenvalues simultaneoulsy become complex across the PT-symmetric 
threshold [18]. In the following, we obtain a-dependent generalizations of these results. 
When a > 0, the bandwidth of a clean lattice scales as A' a (N) ~ N a ; to be consistent with 
the a = case, we use 7/A a as the dimensionless measure of the impurity strength where 
A Q = A' a /4. 

The left panel in Fig. [l] shows the typical phase diagram for an even lattice (N = 100) in 
the (7/A a ,/i) plane as a function of the hopping exponent a > J2JJ. The PT-symmetric 
threshold, below which all eigenvalues are real, decreases with increasing distance between 
the impurities, d = 1 + N(l — 2fi), but approaches a universal critical value 7pt/A q = 1 
when the impurities are the closest, fj, = 1/2. When iV is odd (not shown), the fragile nature 
of the PT-symmetric phase is maintained, but the critical strength for closest impurities is 
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FIG. 1. (color online) a) Typical phase diagram of the Hamiltonian H a as a function of impurity 
strength 7/A Q and fractional impurity position /i = tuq/N for an N = 100 site lattice. The 
region in the (7/A Q ,/u) plane below each curve represents the 'PT-symmetric phase with purely 
real eigenvalues; recall that the quarter-bandwidth A a ~ N a . The T^T-symmetric phase is robust 
at = 1/2. The results are similar for an odd N, except that the critical strength when the 
impurities are closest is jpt = A a /2. b) Scaling of the dimensionless critical impurity strength 
7pt(-/V, /x)/A a for an a = 1 lattice. When the impurities are farthest, jpt ~ 1/vN (blue circles). 
For an intermediate position, [i = 1/4, the critical ^pt ~ iV -1 / 3 shows sub-linear behavior when 
plotted as a function of 1/y/N (red squares). For fi = 1/2, the approach to the asymptotic limit 
jpt{N) — "Ypt{N — > oo) ~ 1/N, plotted here, shows a quadratic behavior as a function of 1/VN 
(beige diamonds). Thus, the "fragile" phase of an a > lattice is significantly stronger than its 
a = counterpart. 

reduced by a factor of two, jpt = A Q /2. In both cases, when a > 0, the critical impurity 
strength 7p-r(A i ) for a generic /i is appreciable compared to its respective universal value at 
// ~ 1/2. Thus, the robust nature of the PT-symmetric phase at fi = 1/N and its extremely 
fragile nature at generic values of fi (a = 0, blue circles) are anomalous features of a lattice 
with constant hopping [18J. 

The right panel in Fig. [l] shows typical scaling of the critical impurity strength jpt^N, fx) 
for /i = {1/N, 1/4, 1/2} in a lattice with a — 1. When the impurities are farthest, /i = 1/N, 
the dimensionless critical impurity strength vanishes as 'Jpt(N) ~ N~ l l 2 (blue circles); when 
jj = 1/4, it vanishes as ^pt{N) ~ iV -1 / 3 (red squares); and when /i = 1/2, it approaches the 
nonzero critical value as 7p T (iV) — ^pt^N — >■ oo) ~ iV -1 (beige diamonds). Note that these 



results are in stark contrast with the a = case where 7pt(^0 ~ N' 1 for all generic /x [18]. 
Thus, although the PT-symmetric phase of a non-uniform lattice for a generic impurity 
position is, in principle, algebraically fragile, the critical impurity strength is still appreciably 
large for experimentally relevant numbers (of waveguides); 7p;r(A r )/A Q , ~ 0.5 — 0.3 for N < 
100. 

Next, we consider the emergence of complex eigenvalues [20]. For a tight-binding Hamil- 
tonian with a purely imaginary potential, if E is an eigenvalue, then so are — E and E* [22J. 
For a given impurity position itlq, the first (four) complex eigenvalues emerge as two adja- 
cent energy levels (-E, —E + S) become degenerate, and so do their positive counterparts. 
Thus, as the fractional impurity position \x changes from the edge of the lattice (// = 1/N) 
to its center (/i « 1/2), the position of the first complex eigenvalue changes from the center 
of the energy band, E = 0, to the its edges, E ~ ±2A a . For a given /i, when the impurity 
strength 7/A Q increases, we find that the maximum number of complex eigenvalues is still 
given by 2/aN = 2mo < N, and that the process is sequential except in an even lattice at 

= 1/2 M- 

Figure [2] shows this generic even-odd effect for the case of closest impurities in a lattice 
with a = 2. For iV = 20 (blue squares), when 7/A a = 1.08, just above the even-lattice 
PT-symmetry breaking threshold, we see that all N eigenvalues simultaneously become 
complex; on the other hand, for N = 21 (red circles), when 7/A a = 0.63, just above the 
odd-lattice PT-symmetry breaking threshold, four eigenvalues near the band-edges become 
complex and the rest remain real. This difference has dramatic consequences. 
Wavepacket Evolution: We now explore the signatures of sequential, or simultaneous, VT 
symmetry breaking in the time-evolution of a single-particle wavepacket. It is straightfor- 
ward to obtain the site- and time-dependent intensity \(k\ijj(t)\ 2 where \ip(t)) = G(t)\ip(0)), 
G(t) pq = [exp(—iH a t/H)] pq is the time-evolution operator, h = h/(2ir) is the (scaled) Planck's 
constant, and (k\ifj(0)) denotes the initial state. We remind the reader that since the Hamil- 
tonian H a is not Hermitian, G(t) is not a unitary matrix. With an open (quantum) system in 
mind, we use the standard inner-product to obtain the time-dependent intensity; therefore, 
the maximum intensity at a site, for a normalized initial state, can exceed unity [23] . 

Figure [3] shows the evolution of the site- and time-dependent intensity across the VT- 
symmetric phase boundary for a particle initially localized on the first site, (^^(O)) = 5^1, in 
an a = 2 lattice. The vertical axis in each panel represents the site index. The top two panels 
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FIG. 2. (color online) Degrees of "P7~-symmetry breaking for an a = 2 lattice with closest impurities 
in the odd (N = 21, red circles) and even (TV = 20, blue squares) cases. The complex eigenvalues are 
scaled by half-bandwidth 2A a , and are obtained for impurity strengths just above their respective 
'PT-symmetry breaking thresholds: 7/A Q = 0.63 for N = 21 and 7/A Q = 1.08 for N = 20. When 
N is odd (red circles), we see that four eigenvalues, near the top and the bottom of the energy 
band, become complex while the remaining (N — 4) remain real; when TV is even (blue squares), 
all ./V eigenvalues simultaneously develop finite imaginary parts. This contrast becomes large for 
iV» 1. 



show the predicted intensity profile for an odd lattice (N = 21) with closest PT-symmetric 
impurities ±27, at positions (10, 12). The first panel shows intensity oscillations that occur 
as the particle bounces from one end of the lattice to the other while encountering the 
(gain and loss) impurities with strength 7/A a = 0.626 at the center. When 7/A Q = 0.627, 
the PT-symmetry is broken with four complex eigenvalues (Fig. [2|. We see that over the 
time-scale shown in Fig. [3j the intensity distribution does not change significantly. The 
bottom two panels are for an N = 20 lattice with impurities at positions (10,11). The 
third panel shows the intensity oscillations that occur in the "PT-symmetric phase when 
7/A a = 1.070. When the impurity strength exceeds the critical value, 7/A a = 0.174, all 
N eigenvalues become complex, thus maximally breaking the VT symmetry (Fig. [2]). This 
property is manifest in the predicted intensity profile in the fourth, bottom-most panel. 
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FIG. 3. (color online) Site- and time-dependent intensity profile for a particle initially localized 
on the first site in an a = 2 lattice. The vertical axis in each panel represents the site index. 
The top two panels are for an odd lattice, N = 21, with impurities ±ij at positions (10, 12). 
They show only a small change in the intensity profile as 7/A Q is increased, from 0.626 to 0.627, 
across the 'PT-symmetry breaking threshold and four eigenvalues become complex. The bottom 
two panels are for an even lattice, ./V = 20, with nearest neighbor impurities at positions (10, 11). 
They show a dramatic change in the intensity profile that occurs when 7/A Q is increased, from 
1.070 to 1.074, across the PT-symmetry breaking threshold and all N eigenvalues simultaneously 
become complex. This even-odd effect for /x ~ 1/2 is independent of N 3> 1 or a. 

Note that we have scaled the color-bars in the bottom two panels to an identical range, to 
simplify comparison while maintaining the visibility of structures in the intensity profile; the 
maximum intensity in the bottom- most panel is ~ 10 4 . 

Lastly, motivated by the recent experimental setup [T7j, we consider time evolution in a 
system where the impurities are given by (7 e (£), 7l) at positions (m , rh ) respectively, where 
7l > is the loss and 7 e (£) = 7g(^) — 1l is the effective, time-dependent, gain [IT] . We choose 
7g(^) = 27l(1 — e~*/ T ) and 7^ = jpt + + , so that when t/r ^> 1, the system is (barely) in 
the broken PT-symmetric phase. Figure [4] shows the remarkable even-odd effect that occurs 
due to different degrees of PT-symmetry breaking. The left (right) panel corresponds to 
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FIG. 4. (color online) Site- and time-dependent intensity profile for a state initially localized on 
the time-dependent-gain impurity, mo = 10, shown by a white dot, in an a = 1 lattice. The 
vertical axis denotes time normalized in units of T a = 1/A a , and the gain-impurity time-constant 
is r/T a = 5. The intensity range for the even lattice (left panel) is larger than that for the odd 
lattice (right panel). The a = 1 lattice spectrum has equidistant eigenvalues \19\ 123], which gives 
rise to periodic behavior seen in both, even and odd, cases at times t/T a < 50. At large times, 
the 'PT-symmetry is maximally broken for the even, N = 20, lattice, leading to commensurate 
real parts for the complex eigenvalues, and persistence of a "periodic" structure in the intensity; in 
contrast, for the odd, N = 21, lattice, the commensuration partially is lost due to the emergence 
of only four complex eigenvalues. The even-odd effect does not depend upon the commensurate 
nature of the lattice spectrum and is valid for any a > 0, including a lattice with constant hopping, 
a = 0. 

results for an a = 1 lattice with N = 20,7 L /A Q = 1.06 (N = 2l,^ L /A a = 0.60), and 
gain-impurity position mo = 10. The initial state (k\tp(0)) = Sk mo , localized at the gain 
site, is indicated by the white dot in both panels. The time-constant r = 5T a = 5/A a is 
chosen so that the PT-symmetry condition, 7 e (t) = 7l, is quickly achieved. We see that 
the overall intensity for the even lattice (left panel) is higher than that for the odd lattice 
(right panel), consistent with a higher number of complex eigenvalues. We emphasize that if 
the impurities are not the closest, the qualitative difference between time-evolutions in even 
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and odd lattices is small; in both cases, then, the PT-symmetry breaks sequentially and 
the initial number of complex eigenvalues is the same. These results are valid for a lattice 
with any a, including a constant hopping, and a generic initial state. They predict that the 
degree of PT -symmetry breaking will be clearly visible in intensity measurements, and that 
even and odd lattices will have dramatically different signatures of a broken PT-symmetry 
when iV > I. 

Discussion: In this paper, we have explored the robustness of the PT-symmetric phase, 
as well as the degree and signatures of its breaking in a lattice with a position-dependent, 
parity symmetric, hopping function t a and a pair of imaginary, PT-symmetric impurities. 
We have shown that when the hopping is maximum near the center of the lattice, a > 0, 
the PT-symmetric phase lasts over an appreciable range of impurity strength, and that it is 
robust when the impurities are closest. The nature of PT-symmetry breaking transition at 
the robust point in even and odd lattices is different: when iV ^> 1 is even, all eigenvalues 
simultaneously become complex, whereas when N ^> 1 is odd, only four of them initially 
do. We have predicted that this difference leads to clear signatures in the intensity profile 
of a single-particle wavepacket. 

In these calculations, we have ignored on-site disorder effects that, in coupled optical 
waveguides, arise due to variation in the real part of the refractive index. The effect of 
disorder on the critical impurity strength, as well as the interplay between the disorder, that 
tends to localize the particle to its initial position, and an imaginary potential, that tends to 
localize the particle at the impurity location, will deepen our understanding of these lattice 
models. 

Our results suggest that PT-symmetry breaking in a (non-uniform) lattice with N ^> 1 
sites will be accompanied by remarkable phenomena, such as the even-odd effect predicted 
here, with no counterparts in the N = 2 system [T6| ITT] . 
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